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Abstract. The simple leaky-box model of propaga- 
tion of cosmic rays in the Galaxy is quite suitable for 
handling data on cosmic ray nuclei energy spectra 
and composition at E » 1 GeV [1,2]. In the leaky- 
box model a full information about cosmic ray prop- 
agation in Galaxy is compressed to the single param- 
eter - escape length, X e , characterizing mean grams 
of a matter passed by cosmic rays from sources to 
the Earth. In this paper we analyze the world data 
on proton and iron cosmic ray spectra collected in 
the past (HEAO, CRN et al.) and in series of recent 
electronic experiments (ATIC, CREAM, AMS, BESS, 
Tracer et al.) and obtain the rigidity dependence of 
escape length, X e (R) =~ # 47± ~ 03 , from the 
measured rigidity dependence of the protons/iron 
ratio. It quite agrees with the one estimated in 
standard manner from the secondary/primary nuclei 
ratio. But at R > 300 GV the behavior of X e (R) 
distinctly changes, that can (variant of explanation) 
point out to the change of proton/iron ratio in cosmic 
ray sources. 
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I. Introduction 

In spite of apparent absence of physical background 
for the leaky-box model, where the transport of energetic 
particles is described by introducing the mean escape 
time of cosmic rays from the Galaxy, and the cosmic- 
ray density, the source density, the gas density do not 
depend on coordinates, it may be applied to the study of 
diffusion and nuclear spallation of stable nuclei attended 
by the production of secondary relativistic nuclei in the 
interstellar gas JT] (2). This can be explained by the 
concentration of cosmic rays sources and the interstellar 
gas in a relatively thin galactic disk immersed in the flat 
but thick cosmic-ray halo |T|. The spatial distribution 
in the low-density halo is the same for different stable 
nuclei because of the negligible nuclear spallation. The 
calculated intensities of stable nuclei for an observer 
at the galactic disk look as corresponding leaky-box 
expressions even for nuclei with large cross sections and 
all propagation is described by some escape length of 



cosmic rays from the Galaxy X e (measured in g/cm 2 ) 
that is a function of particle rigidity R. This important 
parameter means grams of a matter passed by cosmic 
rays from sources to the Earth in average. In [2| there 
was proposed a way to estimate the analogous parameter 
in the diffusion model - effective values of X e f. It was 
shown that leaky-box is a good approximation to the 
widely known and used for the interpretation of various 
cosmic-ray data basic GALPROP (Galactic Propagation) 
model [3 1 with 1% accuracy for all nuclei. In [2| values 
of X e f were found for the three set of GALPROP 
parameters, corresponding to the three models of the 
propagations. 1) Plain diffusion model (PD) 

X ef = 19/3 3 g/cm 2 R < 3 GV, 

X ef = 19/3 3 {R/3 GV)-°- 6 g/cm 2 , R > 3 GV; (1) 

2) Diffision with acceleration (DR) 

X ef = 7.2{R/ 3GV)~°' 34 g/cm 2 , R > 30 GV; (2) 

3) Diffusion with damping (DRD) 

X ef = 13(i?/3 GV)-°- 5 g/cm 2 , R > 10 GV (3) 

For the comparison we present also the widely used 
approximation (4): 

X ef = 11.8/3(i?/4.9 GVr - 54 g/cm 2 ,Fi- > 4.9 GV 

(4) 

The leaky-box approach due to its simplicity until 
very recent was used by many authors. 

Intensity I a of CR nuclei with mass number A near 
the Earth is connected with their spectrum in a source 
Qyi(if neglecting the energy losses at small energies 
and the contribution from the fragments produced from 
heavier nuclei) by the simple equation: 



Ia 



4irp 



l 



(5) 
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TABLE I: The fit of the ratio of Ijj/Ip e — fit by the polynomial function of the fifth order. Averaged points of 
the ratio Ijj jJ-Fe ~ points are calculated in 14 R-bins with the corresponding errors 



R, GV 


3.16 


5 


8 


12.6 


20 


31.6 


50 


89 


178 


355 


708 


1412 


2818 


Ip/Ife - fit 


2907 


2313 


1875 


1564 


1347 


1197 


1095 


1012 


952 


914 


889 


877 


894 


Errors 


249 


212 


105 


47 


43 


101 


41 


56 


65 


84 


143 


166 


150 


Ip/Ife — points 


2432 


2281 


1748 


1676 


1235 


1203 


1067 


991 


842 


1053 


955 


976 


1274 


Errors 


249 


212 


105 


50 


40 


101 


40 


54 


57 


97 


154 


185 


447 



where p - gas density, - interaction length for the 
nucleus with mass number A. (calculated from GAL- 
PROP). For protons X? ~ 110 g/cm 2 at 1 GeV and it 
decreases to 70 g/cm 2 at 10000 GV, that is X ef < 
from whence the well known equation arises: 

Ih ~ Qh * Xef, lobs — Isour + 

where a is the slope of X e f dependence, j b s , J SO ur - 
the slopes of observed and source spectra correspond- 
ingly in the case of power-like laws. 

But while interaction length for the iron nuclei is 
much smaller X.^f ~ 2.7 g/cm 2 , there is the rigidity 
region Rmin + Rmax for Fe nuclei where 

iFe ~ QFe * ^in, lobs ~ Isour, 

reaching the asymptotic value (proton and Fe spectra are 
parallel), 7 o6s = j sour + a when X ef at R max becomes 
much smaller than Xf^ . 

We see that one can try to estimate the rigidity depen- 
dence of X e f from the ratio of any two spectra in the 
interval R m i n -j- R m ax, but with essential reservations: 

- to bear in mind that the intrinsic property of the leaky- 
box model is the independence of X e on cross section 
and as a result on the type of nucleus; 

- to assume that all types of nuclei are produced in the 
same type of sources and the chemical composition of 
accelerated particles in the sources does not depend on 
rigidity in the investigated interval - 

Q Fe (R)/Q H (R) = const; 

- to use nuclei with significantly different interaction 
lengths, besides, the contribution of fragments among 
them should be negligible; 

- to find way to estimate Qp e (R)/QH(R)- 

Proton and iron spectra are most suitable ones for this 
task. 

The main question is: if available in our days experi- 
mental data are enough for this task? 

II. Experimental data 

As we are going to analyze spectra at fixed rigidity, 
we need the spectral data measured at remote intervals 
of energy per particle. Interval pointed out in title, 3- 
3000 GV, corresponds to the energy per particle interval 
2.3 GeV - 3 TeV for protons and 44 GeV - 82 TeV 
for iron nuclei. So we are doomed to use data from 
different experiments. Moreover only with appearance 
of new data on iron nuclei from the experiment Tracer 
|[T9l . measured in uniquely wide energy range (E part =3Q 



GeV-80 TeV), this task probably can be solved. Data 
obtained in ATIC2 experiment ifTTll are very important 
also, because they fill up the gap between data ob- 
tained by magnetic spectrometers and ones obtained with 
calorimeters and emulsions. 

We include in the consideration experimental points, 
which satisfy the conditions: 

a) energy interval lies between 3 and 10000 GV, 

b) errors don't exceed 30%, 

c) calibration on accelerators is done (that is why Sokol 
experiment was excluded). 

Fig. 1 and Fig. 2 represent proton and Fe spectra 
measured in different experiments. The points satisfy- 
ing above enumerated conditions are denoted by filled 
symbols. 

Proton data are got from the following experiments: 
MASS91 0, BESS-98 [6), 0BESSO2, AMS ©, 0, 
IflOl . ATIC02 flU, CREAM [12], JACEE [13], Ifl4l, 
RUNJOB |Q3). Fe nuclei are from: HEAO-3 JTgl, CRN 
ifTTl . TRACER (Til, ED, ED, Sanriku J20). Data were 
reduced to fixed modulation potential <f> = 800 MV by 
calculation similar to iflOl . 

Fortunately points satisfying above mentioned condi- 
tions are little scattered. For further work and conve- 
nience one needs to approximated these dependencies 
by any way. In the Fig. 1 and 2 we denote by thick 
line the best fit of corresponding scatter points, that is 
a polynomial function of the fifth order. In the Table 1 
this fit is denoted as " Ijj / If& — fit " with errors in the 
upper two rows. 

Moreover, we calculated averaged points in the every 
from 14 R-bins getting errors which include statistical 
significance and a divergence of experimental points 
fallen in the fixed bin (thin lines with errors in Fig. 
1,2.). In the Table 1 averaged points are denoted as " 
Ih/Ife — points " and they are presented in the lower 
two rows. 

III. Results and discussion 

The leaky-box Ed. (5) written for the proton and 
Fe spectra allows us to express the value of effective 
escape length X e f by means of the ratio of proton/iron 
measured spectra Ih/Fc = Ih / Ifc{R) (from Table 1), 
the source chemical composition Qh/Fc = Qh/Qfc 
(which does not depend on rigidity and should be 
determined any how) and known values of Xf n e , X^: 

Y lH/Fe(R)/QH/Fe - 1 , , 
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Fig. 1: Proton spectra measured in different experiments (for ref. see text), reduced to the modulation potential 
$ = 800 MV; open symbols - points with more than 30%; thick line - the 5th order polynomial fit, thin lines - 
averaged points 



LO 

£ 10 2 

> 
CD 



10' 



CO 



CM 



^ 10° 

CM 

T3 




J 



• TRACER B HEAO-3 

★ , * CRN —a— A Sanriku 

▼ ATIC02 □ JACEE, O RUNJOB, 
Fit, — ° — averaged points 



T 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I I 



10 1 



1 1 1 — I — I I I I I I 

10 2 10 3 
R, GV 



10" 



Fig. 2: Spectra of iron nuclei; open symbols - points with more than 30%; thick line - the 5th order polynomial 
fit, thin lines - averaged points 



The Eq. (6) will work obviously only in the range 
of rigidities R m in — Rmax where X e f is not much 
smaller than Xf r f: at X e f — > I H / Fe (R) -> Qu/Fe- 
To demonstrate expected application range of Eq. (6) 
we choose the critical value R max corresponding to 
X ef = 0.1Xf n e ~ 0.27g/cm 2 . It means that the region 
X e f < 0.27g/cm 2 and R > R ma x is the asymptotic 
region where spectra become parallel and not sensitive to 
the X e f. For three models embedded in GALPROP (see 
Introduction) the corresponding values of R m ax look 
like below: 

1) R max = 3.6 TV for PD model(a = -0.6), 

2) R max = 6.9 TV for DRD model(a = -0.5), 



3) R max = 46.8 TV for DR model(a = -0.34), 

4) R max = 69000 TV for RD model(a = -0.2) for the 
very flat energy dependence X e f — 10R~ 0/2 considered 
in Ell being chosen for the explanation of the "knee" 
in PCR by the change of propagation mechanism. 

So Rmax comparable with the experimental maximal 
R = 2.8 TV (see Table 1) may be only for very steep 
X e f dependencies, as in the cases of PD or DRD models. 
In this case we determine the value of Qn/Fe from the 
asymptotic ratio of In/FeiR) taking into account minor 
corrections. We estimate Qn/Fe — 800 — 1000. 

The variant of calculation by the Eq. (6) for the 
experimental values of Ijx/p e {R) from the Table 1 and 
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Qff/Fe = 800 give the next approximation: X e f(R) ~ 
4.6 * (R/5 GV)~ - 65 , that is much lower than the 
expected dependence X e f estimated by the B/C ratio 
(l)-(4). This means that our assumption of steep X e f 
decrease with energy R> 3 TV (as R~ Q5 ~- 6 ) is not 
fully correct. 

For the analysis of the set of flat dependencies 
X e f(R) where R ma x 3> 3 TV and there is no possi- 
bility to estimate the asymptotic values I H / Fe (R) it is 
proposed to "normalize" X e f at R — 50 GV by the one 
obtained from B/C measurements. The values of X^j/ C 
obtained from B/C - ratio occur between 2.8 g/cm 2 (DR 
model - (2)) and 3.45 g/cm 2 (PD model (1)). The point 
R = 50 GV was chosen because the large body of data 
on B/C are in a good coincidence in this point 11221 . 
from the other side at the R = 50 GV the contribution 
of reacceleration processes surely should be small. 

Then Qh/Fc value could be determined 
using X^ /C (50GV) « 2.8 -3.45 g/cm 2 , and 
I H /Fe(R){^GV)= =1067+40 (from Tabl. 1). By 
Eq. (6) we get the corresponding interval for Qn/Fe = 
490 - 535. Substituting obtained limit values of Qn/Fe 
and measured regidity dependence of proton/iron 
spectra from Tabl. 1 in Eq. (6) we calculate two 
variants of rigidity dependencies of X e f(R), presented 
in Fig. 3. The upper one (with Qn/Fe =490) could be 
approximated in the interval R = 3 — 300 GV as 
X ef {R) ~ 3.46 * (R/5 GV)-°- 47±a03 .. The lower one 
(with Qn/Fe =535) could be approximated as 
X ef (R) ~ 2.7 * (R/5 GV)-°- 50±0 - 03 . 
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Fig. 3: X e f(R) dependencies calculated from the mea- 
sured proton/iron spectra (lowest two rows in Tab. 1) 
for the Q H / Fe =535 (black circles), Qh/_fe = 490 (open 
circles). Thin lines - the approximations by power law. 

Up to the R ~ 300 GV the obtained from protons/iron 
ratio rigidity dependence of escape length X e f(R) can 
be approximated by the power-like law: X e f(R) ~ 
i?-o.47±-o.03^ that is in a good agreement with the one 
estimated in standard manner from secondary/primary 
nuclei ratio (l)-(4). But at R > 300 GV the behavior of 
X e f(R) changes significantly. It is caused by the "im- 
proper" behavior of the slopes of proton and iron spectra 
at R > 100 GV, 7Fe = 2.71±0.03 [19| and 7^=2.63 ± 
0.03 ifTTIl . They are in a striking contradiction with the 



expected values: the Fe spectrum should be flatter than 
proton spectrum (see Eq.(5)). 

Here it is worth noting that in the region 300 GV-3 TV 
(Fig.l, 2) the principal contribution for protons comes 
from ATIC2 data [ 1 1 1 and for iron nuclei it comes from 
the Tracer data 1 19 1. 

So if we believe that both these important experiments 
are reliable, we should inevitably conclude that in the 
range R > 300 GV there is a change of proton/iron 
ratio in cosmic ray sources. Some authors of ATIC have 
come to this conclusion already 11231 . 

In summary we would like to stress that the region 
300 GV - 10 TV continues to be of importance, and the 
result that spectrum of iron nuclei is flatter than proton 
spectrum needs to be confirmed in other experiments. 

The work is supported by RFBR grant 10-02- 
01443-a, V.S.Ptuskin - RFBR grant 10-02-001 10a. 
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